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Abstract. We study the statistics of charge transport in a chaotic cavity attached to 
external reservoirs by two openings of different size which transmit non-equal number 
of quantum channels. An exact formula for the cumulant generating function has 
been derived by means of the Keldysh-Green function technique within the circuit 
theory of mesoscopic transport. The derived formula determines the full counting 
statistics of charge transport, i.e., the probability distribution and all-order cumulants 
of current noise. It is found that, for asymmetric cavities, in contrast to other 
mesoscopic systems, the third-order cumulant changes the sign at high biases. This 
effect is attributed to the skewness of the distribution of transmission eigenvalues with 
respect to forward/backward scattering. For a symmetric cavity we find that the third 
cumulant approaches a voltage-independent constant proportional to the temperature 
and the number of quantum channels in the leads. 
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1. Introduction 

Statistical properties of charge transfer in mesoscopic conductors have recently received 
rapidly growing attention. A powerful theoretical approach to the fluctuation problem 
is based on the concept of full counting statistics (FCS) [1 . It gives the probability 
to transmit a given number of charges from one terminal of a mesoscopic conductor to 
another terminal during the observation time. From the knowledge of these probabilities 
one can readily derive not only average current and noise, but all higher cumulants of 
current fluctuations and extreme-value statistics. Thus FCS represents the complete 
information which can be gained from charge counting in a transport process. 

Since the introduction of FCS for electron counting in mesoscopic conductors 
by Levitov and co-workers jT], the approach has been successfully applied to various 
systems including quantum point contacts, diffusive conductors, cavities, metal and 
superconductor junctions 2j. Until very recently, the experimental measurements 
were focused mainly on shot noise — the second cumulant of the fluctuating current 
[3]. Experimental investigation of FCS has been opened up by measurements of the 
third cumulant of current fluctuations jU IH]. Since the third cumulant (odd-order 
cumulant) changes sign under a time reversal transformation, it is not masked by 
thermal fluctuations and can be observed at relatively high temperatures [H]. In contrast, 
accurate measurements of the second cumulant or shot noise (even-order cumulant) is 
difficult because of the requirement to maintain a low temperature at a high applied 
voltage. 

In this paper we address the FCS of an asymmetric chaotic cavity with two 
attached leads supporting non-equal number of quantum channels. Previous studies 
on asymmetric chaotic cavities include the following. The second cumulant has been 
calculated by Nazarov [7j by means of a circuit theory, by Beenakker [S] using random- 
matrix theory (RMT), by Blanter and Sukhorukov semiclassically [Oj, and by Schanz et 
al [ini using quantum graphs. The third and fourth cumulants of current were obtained 
by Blanter et al jlT]. The authors found a discrepancy between semiclassical and 
RMT results for high-order cumulants starting from the fourth cumulant. Subsequently, 
Nagaev et al were able to explain the discrepancy by taking into account "cascade 
correlations" ^3] between the cumulants within a semiclassical approach. It was pointed 
out, that evaluation of high-order cumulants goes beyond the validity of the standard 
Boltzmann-Langevin scheme: one has to use cascade diagrams, or to introduce non- 
local noise correlators in the Boltzmann-Langevin equation It should be noted 
that the results on the third and forth cumulants of current fluctuations |^ were 
obtained in the zero temperature limit only. Experimental measurements of the second 
cumulant (shot-noise) on chaotic cavities were performed by Oberholzer et al I16j . 
Recently, a stochastic path integral approach has been formulated by Pilgram et al 
PT] and Jordan et al [18 , which is efficient to study FCS for arbitrary mesoscopic 
semiclassical systems including stochastic networks. The authors ^Hj have derived the 
diagrammatic rules to calculate the high-order cumulants of current fluctuations which 
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take into account the Nagaev's cascade correction terms ^21- This approach was apphed 
to calculate the cumulants of charge fluctuations in asymmetric cavities in the regime 
of strong electron-electron scattering [TH Hn] . Pilgram and Biittiker studied the effect 
of screening on fluctuations of charge inside the cavity 

In this paper, we use the extended Keldysh-Green function technique within the 
circuit theory introduced by Nazarov |2I]. By using this technique, we have obtained the 
cumulant generating function (CGF) in exact analytical form which contains complete 
information about all-order cumulants of current fluctuations in asymmetric chaotic 
cavities at arbitrary bias-to-temperature ratio. The analytical formula for CGF is 
expressed in terms of the anticommutator of the matrix Green functions in the leads. 
In this form it can be applied to calculate all-order cumulants of current noise and 
extreme- value statistics for different types of the contacts. In section |2l we present the 
derivation of CGF for asymmetric cavities and show the correspondence with the results 
obtained earlier in the symmetric limit. The distribution function of charge transfer is 
obtained in section IHl The flrst three cumulants of current fluctuations are analyzed in 
section m for arbitrary bias-to-temperature ratio. The third-order cumulant was found 
to be linear in voltage at both low and high biases, but, in contrast to other mesoscopic 
systems, it changes the sign at high biases for any asymmetry of the cavity except the 
ideally symmetric case. For the symmetric cavity we found a saturation of the third 
cumulant with bias at the value proportional to the temperature and the number of 
channels in the leads. We explain the change of sign for the third cumulant as a feature 
of skewness of the distribution function of transmission eigenvalues. This is conflrmed by 
the analysis of the temperature dependent third cumulant for an arbitrary mesoscopic 
conductor in sectional In the Appendix we show the correspondence between the CGF 
obtained in the present paper and the RMT results in the zero temperature limit. 

2. Cumulant generating function 

A chaotic cavity is a conducting island of irregular shape connected to external leads 
by two quantum point contacts (flgure HJ. The openings are assumed to be small 
compared to the size of the cavity. The leads are connected to electron reservoirs kept 
in local thermal equilibrium with chemical potentials fi^ and fin. The difference in 
chemical potentials f^L — fJ^R = gives rise to a current / flowing through the cavity. 
For simplicity, the quantum point contacts are assumed to be completely transparent, 
i.e., their widths are adjusted to conductance plateaus [IH]. We disregard possible 
Coulomb blockade effects, which is justifled if the conductances of the leads Gl = GqNl, 
Gr = GoNn are assumed to be much larger than the conductance quantum Go = e^/nh. 
This means a multichannel transport regime and negligible weak localization effects. 
The inelastic relaxation inside the cavity is disregarded, so the energy of each electron 
entering the cavity is conserved. Since we are interested only in low-frequency limit 
of current fluctuations, the screening effects in the cavity are neglected. It is also 
assumed that the measurement circuit has negligible impedance in order to avoid the 
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environmental effects on the statistics of current fluctuations 

We apply a circuit theory of FCS based on equations for the Keldysh-Green 
functions [2Ij. Its main idea is to represent a mesoscopic system by a set of discrete 
elements: terminals, nodes and connectors. In contrast to conventional classical circuits, 
in which the nodes are described by a scalar electric potential, here the nodes are 
characterized by matrix Green functions. The matrix currents flowing between the 
nodes and obeying matrix analogs of Kirchoff rules contain all current correlators 

[211 ESI El Ea EE]. 

The circuit theory representation of a chaotic cavity is depicted in figure [T] (right 
panel). It consists of two terminals (the left and right reservoirs), a central node (the 
cavity) and two connectors (the left and right point contacts). The terminals and the 
node are described by 2x2 matrix Green functions in Keldysh space: Gl, Gr, and Gc, 
respectively. The terminal matrices are assumed to be known and they are given by 

(5^ = e^^-/2G0e-'^-/2, Gr = G% (2.1) 

where the counting rotation is applied to the left terminal, is the third Pauli matrix 
and the energy dependent Fermi function /^(e) = {1 + exp[(£: — fik)/kTk]}^^ accounts 
for the electrochemical potentials /i^ and the temperatures in the reservoirs. For a 
two-terminal conductor, it is sufficient to introduce only one counting field A by gauge 
transformation in the left terminal. Since we focus only to low-frequency limit of current 
fluctuations, A is taken as a time independent constant. 

The matrix current between each connector k and a central node is given by [SB] 

J = V Tn [G'fc, Gc] ,r) n\ 

' t-4 + T„^({Gfc,GJ-2)- 

Each connector k is characterized by a set of transmission eigenvalues associated 
with n quantum channels at energy e. For fully transparent ballistic point contacts 
= 1 for n < Nk and = for n > A^^, that corresponds to having Nl and Nr open 
channels. In general A'^^: 7^ Nr. Thus the matrix current conservation 1^ + Ir = gives 

N, [^"J ^ + Nr ^^^^J ^ = 0. (2.4) 
^2 + {Gi,GJ ''2 + {Gr,Gc} ^ ' 

Here, the only unknown function is the matrix Gc-, which is isotropic within the cavity 
under a quantum chaotic regime. Since the terminal matrices Gk and the node matrix 
Gc are traceless, they can be expanded over the Pauli matrices |25j: Gk = VkT, 
Gc = VcT, where Vk and Vc are three-dimensional vectors, and r = {fx,fy,Tz) is the 
vector of the Pauli matrices. With the use of this parametrization the anticommutators 
{Gk,Gc} = 2vk ■ Vc and {Gl,Gr} = 2vi ■ vr are scalar functions. Introducing the 
notation Gy, = VlGl + PrGr. with 
_ Nk 
2 + {Gfc,GJ' 



Full counting statistics of a chaotic cavity 



5 



we rewrite equation (j2.4p as [Gj],Gc] = 0, from which we obtain PlVl + PrVr = cVc, 
where c is the constant to be found. Using the normahzation conditions v^=v^=l and 
denoting = (A^'^, + Nfi)/2, we obtain the set of equations 

c'=pI+pI + PlPr {Gl, Gr}, (2.6) 

2pL = iV+ - cNn/N+, (2.7) 

2pn = N+~cNL/N^, (2.8) 
from which the constant c can be found exphcitly: 

c = N+{1 + VT^)~\ (2.9) 



{Nl + NrY \ 2 + {Gl,Gr}I' 
Our aim is to find the CGF S{X) related to the matrix current by 



(2.10) 



-^uj-MU (2.11) 

where t is the observation time, which is assumed to be much longer than both the 
correlation time of current fluctuations and the characteristic time e/I with I the time- 
averaged current. The latter condition ensures that the number of electrons transfered 
during the time t is much greater than 1. Substituting II = Pl [Gl,Gc] and doing 
matrix transformations we obtain 

with 

A=^-^{GL,Gn}. (2.13) 

To find the CGF from equation (j2.12p explicitly, we have to integrate A over A. To this 
end all the variables in equation ()2.13|) can be expressed in terms of the anticommutator 
g = {Gl, Gr} by means of equations (jTT)|) - (jZTn|l . We get 



1 — \/l — a dQ 

where the function a(A) is determined by equation (j2.1(jp . Then we use the parametriza- 
tion a = 4x(l — x), that leads to 

(1 — 2a;) dx 

A = 2N+ ^ T^Tn' 2.15 

a;[l — x(l — xji' OA 



where 



NlNr 

" (2.16) 



(NL + NRy 

is the constant which will appear below as the Fano factor. Integrating equation ()2.12j] 
over A by making use of equation ()2.15|) . we obtain the final formula for the CGF 

S(A) = t (iVi + Nr) [ ^ SeiX), (2.17) 
J Tin 

S,{X) = aL In f-^V) + fr^Vl ' (2.18) 

\l-aRlxJ \l-aL/x) 
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with 

x(A) = [1 + v'l - a(A)]/2, (2.19) 

c^k= , k = L,R. (2.20) 

It is important that the analytical result for the CGF is expressed through the 
anticommutator {Gl, Gr} of arbitrary (traceless) 2x2 matrix Green functions at the 
terminals. Therefore the obtained formula is valid for different types of contacts, f 
For normal metals with the help of equation (j2.2|) . the anticommutator Q reads 

g = 2 + 4[h{l - fR){e'' - 1) + fnil - h){e-'^ - 1)]. (2.21) 

It is worth noticing that our results are in agreement with the RMT theory. As we 
show in the Appendix, the distribution of transmission eigenvalues p{T) obtained from 
equation ()2.17|) in the zero temperature limit gives the RMT result. We have also 
verified the correspondence of our results with the stochastic path integral approach. 
The saddle-point solutions of the action in jlTj can be found not only in the symmetric 
limit, but quite generally for an arbitrary asymmetric cavity, and we have checked that 
the result is identical to equation (j2.17|) . 

The semiclassical electron distribution function inside the cavity can be obtained 
from the matrix Green function in the central node 

a = clGl + anGR + [a;(A) - 1]{Gl + Gr). (2.22) 

In the limit A — > we obtain Gc = o^lG^ + aRG% which with the help of equation ()2.2p 
gives the distribution function in the cavity fc = ctLfi + oir/r derived earlier by other 
methods p. 

In the symmetric limit Nl=Nr=N, aL=oiR=l/2, the CGF in equation ()2.17|1 is 
simplified to 



Sf'^iX) = ln[(2 + V2 + g)/% (2.23) 

which corresponds to the result derived by Belzig ^| with the normalization S{0) = 0. 

It is instructive to compare the symmetric and asymmetric cases in the zero 
temperature limit. We present the results for the characteristic function % = e*^: 

/I 1 

X^y^iX) = + -e'V2j ^ (2.24) 

/ s MNl / \ MNu 

,(A) = + ^fp^) U + ^^) . (2.25) 

Here, M = eVt/ nh is the number of attempts per channel, which is the maximal number 
of electrons that can be transferred under the restriction imposed by the exclusion 
principle for fermions ^J. Equation ()2.24|1 has been obtained by Blanter et al from 



f After submission of the first version of this paper to the e- archive [Preprint cond-mat/0403388 
Vanevic and Belzig have proved ^Z] that formula H2.17|l for the CGF is independent of the matrix 
structure and, in particular, it is valid for 4x4 matrices in Keldysh-Nambu space for superconductor 
contacts. 



Full counting statistics of a chaotic cavity 



7 



the RMT theory ^T] and by Pilgram et al by the stochastic path integral method |2(Jj . 
Equation ()2.25|) is an extension to the case of cavities with non-equivalent leads, in which 
the relative numbers of channels Ok are defined by equation ()2.2()j) . The particular form 
of can be attributed to interplay of binomial statistics of the two point contacts 
in series (cavity openings) with inter-channel mixing inside the cavity. Therefore the 
logarithmic functions in equations ()2.17|) and ()2.23|) come from the binomial statistics 
of two quantum point contacts Ij. The function x in equations ()2.17p and ()2.25|) and 
the square root in equation ()2.23|1 account for the inter-channel mixing inside the cavity 

m- 

It is interesting to analyze formula ()2.17|l in the strongly asymmetric limit when 
one opening supports many more channels than the other. Let us suppose N^i < N^. 
Denoting the parameter of asymmetry rj = Nr/Nl -C 1, and expanding the CGF in 
series of r], we get 

S,{\) = (InO V + (r' - 1 - InO + 0(r/3), (2.26) 
e(A) = [2 + ^(A)]/4. 

In the zero temperature limit ^ — e^'*', this simplifies to 

Se{X) = iXr] + (e"'^ -l-iX)rf + 0{Tf). (2.27) 

We observe that the leading-order term ~ is linear in A. This means that the charge 
transport in the leading-order is just the stationary current with no noise. The current 
is determined by the opening with the smallest number of channels I ~ GqVNu. It 
does not mean, however, that there are no current fluctuations at all. The second 
cumulant (noise power) and all higher cumulants are determined by the next order term 
(e~'^ — 1) rf, which depends on both openings. For the noise power, for instance, we 
obtain Pj 2el {Nh/Nl). At finite temperatures the leading-order term ~ rj affects all 
the cumulants of current fluctuations [see equation ()2.26|) ]. 



3. Electron counting statistics 

Having found the CGF S'(A), we obtain a complete information about the transfer of 
charges with time through a chaotic cavity. The full counting statistics is determined by 
the distribution function P{Q) giving the probability of transferring Q charges through 
a cavity during the observation time t. The distribution function is related to the CGF 
by Fourier transform: 

J-TT 277 

This integral can be evaluated in the stationary phase approximation. In the low- 
temperature limit the stationary point is determined by 



{Nl + Nr)M 



iQ. (3.2) 
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The large number of channels Nl^^i ^ 1 and long observation time M ^ 1 provide a 
necessary condition for the observation of a large number of transmitted particles and 
accurate estimate of the integral in the stationary point. By using the CGF ()2.18j) we 
arrive at the pure imaginary value for the stationary point 



The maximal number of electrons that can be transferred through each contact 
separately is given by the product MNk (the number of attempts per channel multiplied 
by the number of channels) with k = L,R . Therefore the charge Q transmitted through 
a cavity is restricted by the limits: < Q < Qmax, where Qmax = MNr and A''^ is the 
number of channels in the narrowest contact. We notice that in this interval of allowed 
values of Q the expression in the logarithm ()3.3|) is strictly positive and the stationary 
point does exist. Evaluation of the integral (jH.lj) in the stationary point (jH.Hj) leads to 
the result 



where Qk = MNk, k = L,R. It is interesting to see that the probability function 



The result of the evaluation of InP(Q) using Eq. ()3.4|) is shown in Fig.Elfor different 
asymmetries of the cavity r] = Nr/Nl. For a symmetric cavity {ri=l) the distribution 
function is symmetric around the average value Q/Qmax = 0.5 and it is in agreement 
with the known result: lnPsym(Q)/Qmax = -2 [ln2 + Qln Q + (1 - Q) ln(l - Q)] with 
Q = Q/Qmax (see, e.g., Ref. [2H1)- When rj is decreased from 1 to 1/100 the cavity 
becomes more and more asymmetric, the maximum P{Q) is shifted from 0.5 to 1 (see 
Fig. 12]). This evolution correponds to the widening of one of the contacts, while the other 
(the narrowest one) is kept fixed. In the strongly asymmetric limit when one opening 
supports many more channels than the other (77 -C 1), the distribution function becomes 
highly asymmetric with a long tail at Q < Q and a sharp drop to zero at Q > Q. The 
maximum of P{Q) approaches the largest value Q — > Qmax determined by the smallest 
opening. This is in agreement with the fact that the conductance of the cavity (the 
first-order cumulant) in this limit depends only on the smallest opening. The width of 
P{Q) determines the second cumulant (the noise power). In the strongly asymmetric 
limit the second cumulant and all higher cumulants are determined by the tail of P{Q) 
at Q < Q described asymptotically as \nP{Q) ~ {Qmax — Q) ^"^{Nr/ Nl) anf therefore 
they depend on both openings. This clearly shows the non-Gaussian character of the 
fluctuations of the transmitted charge in asymmetric cavities. 




(3.3) 



In P(Q) = -2Q In g - ^(g, - Q) \xi{Qk - Q) 



(3.4) 



k 




with 
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4. Cumulants 

The cumulants are found from the CGF by differentiating over the counting field A: 

Ck=[^^ S{X)\x=o. (4.1) 

The first few cumulants characterize the features of the distribution P{Q): the first 
cumulant gives the mean, the second, the dispersion, and the third, the skewness of the 
distribution [Hj. 

Mean current and noise 

Although the first two cumulants for the chaotic cavity are already known, it is 
instructive to see how they emerge in our approach in order to compare with higher 
cumulants. The first two cumulants of current fluctuations determine the mean current 
and noise: 

e 2e2 
/ = ^C7i, Pi = ^C,, (4.2) 

From equations ()2.17|) and ()4.1|) we obtain 
NlNr f de 

I = e——— —Ul- fR), (4.3) 



= / ^ {/.(I - A) + - /.) H- ^^^^(A - h?] ■ (") 

For further calculations we assume the same temperature for all the system, Tl = Tr = 
T, and max(/cT, eV) fik, k = L,R. We get 

/ = GV, (4.5) 

Pi = 4:kTG[l + F (ucothu-l)], (4.6) 

where G = GLGii/{GL + Gr) is the conductance of the cavity, u = eV/{2kT), and F 
is the Fano factor given by equation ()2.16|) . Equation ()4.6|) has also been obtained 
semiclassically jHI EJ. It is seen that the current / is linear in bias for arbitrary 
temperature and therefore contains little information about the transport. In contrast, 
the behavior of the noise power differs in the two limits: 

^ r AkTG[l + {F/l2){eV/kT)\ eV « kT 
^ \ 2eGVF, eV^kT. ^ ' ' 

At low biases (high temperatures) the second cumulant, as an even-order cumulant, 
is dominated by thermal noise 4kTG. To extract information additional to that 
obtained from the mean current, one should go to the shot-noise regime eV ^ kT. 
The measurement of the shot-noise suppression factor Pj/2ql = F, gives in principle 
information about the cavity openings. In these measurements a natural difficulty is 
the requirement to maintain the sample cool under high currents, and to unambiguously 
interpret the linear dependence of Pj over J as a shot noise and not as a thermal noise 
modified by nonlinear conductance jH]. In contrast, the third cumulant is free from this 
difficulty as will be shown in the next section. 
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4-2. Third- order cumulant 

For the third cumulant of current fluctuations 

C/=yC3, (4.8) 

from equation (j2.17|) we find 

+/r(1 - + 272(/i - , (4.9) 
where 71 = 1 — F and 72 = 1 — 2F + 2F^. Integrating over energy, we get 



Ci = ekTG ■ 6F(1 - 2F) 



/^2(1-F) . ^ 

coth u + u \ — — coth^ u 

V3(1-2F) 



(4.10) 



which is valid for arbitrary ratio u = eV/2kT between the bias and temperature. 
Asymptotically, 

^^fe^/F, eV<t:kT 
^ \ eGF[6{l-2F)kT -{l-AF)eV], eV ^ kT. ^' ' 

In the zero temperature limit equation 1)4.111) is in agreement with the results in [mil2j : 
Cj{kT ^ 0) = -e'lF{l - 4F) = -e'l ^^y^T^ . (4.12) 

At high temperatures we observe that close to equilibrium Ci is linear in bias and 
independent of temperature, and hence not masked by the Johnson-Nyquist noise. 
Furthermore, it is proportional to the Fano factor F and therefore the measurements of 
Cj/e^I at arbitrarily small biases, give directly the factor F with no need to subtract the 
thermal noise, in contrast to the case of the second cumulant [compare with equations 
(14.7)) ]. Another important feature is the change of sign of the third cumulant at high 
biases. This behavior is illustrated in figure 121 where the third cumulant is shown as a 
function of bias-to-temperature ratio for different asymmetries of the cavity 77 = Nr/N^. 
[rj is obviously related to the Fano factor, F = r7/(l+?7)^]. For the symmetric case {ri=l), 
Cj approaches the bias-independent constant 

Cf^ieV > kT) = -eOkT = -eGokTN (4.13) 
4 8 

proportional to the temperature T and the number of channels and otherwise 
independent of the parameters of the material. For asymmetric cavities {0 < rj < 1), the 
curves change the sign at eV* /kT = 6(1 — 2F)/(1 — AF) which approaches 6 as r/ ^ 0. 
Equations ()4.7|) and 1)4.111) also demonstrate that with the increase of temperature under 
fixed bias the third cumulant saturates, while the second cumulant diverges with T. 

Next we focus on the voltage dependence of the third cumulant for the asymmetric 
cavity (see figure EJ. At low and high biases, Cj is approximately linear, with positive 
and negative slopes, respectively. These linear dependences can easily be deduced from 
equation ()4.11|) . Note that the slopes, as well as the intersections with the axes, can be 
used in the experiment to determine the parameters of the cavity (the number of channels 
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in each opening). There is a striking similarity of our results with the behavior of Cj 
for a tunnel junction in the electromagnetic environment modelled by two conductors in 
series |22]- For the latter, the temperature dependence and the sign change are caused 
by feedback from the electromagnetic environment. In the case of a cavity, the sign 
change is solely due to internal dynamics. This similarity is probably related to the fact 
that the cavity can also be represented by two conductors (point contacts) in series with 
inter-mode mixing in the node. 

An important point is whether the temperature-dependent third cumulant obtained 
from the circuit theory [equation ()4.9p ] can also be obtained within the semiclassical 
Boltzmann-Langevin picture. The analysis of equations in the work by Nagaev et al 
\12\ shows that the same formula for Ci can only be obtained if one takes into account 
the "cascade corrections". [The authors present the results for Ci at zero temperature 
only.] At zero temperature these corrections are only significant for cumulants starting 
from the forth cumulant, but at finite temperature they are essential for the third 
cumulant too. This means that the third cumulant we have obtained from the circuit 
theory cannot be obtained from the standard Boltzmann-Langevin scheme; one has 
to add the cascade correlations between the cumulants I13j. In addition, it can be 
shown that the bias independent saturation value of C/ for a symmetric cavity [equation 
()4.13|) ] originates completely from the "cascade correction" term. The contribution of 
the minimal correlation term to this value is zero. Thus the experimental measurements 
of the saturation value of the third cumulant would give the evidence of the cascade 
correlations in chaotic cavities and the non-locality of the noise correlators. 

5. Third-order cumulant for an arbitrary mesoscopic conductor 

Having found the change of sign in the third-order cumulant for chaotic cavities, it is 
interesting to ask: (i) is this behavior typical, and (ii) what are the necessary conditions 
to observe the same phenomenon in other types of mesoscopic conductors? Here we 
focus on the effects caused only by internal dynamics due to elastic scattering; the 
environmental effects are neglected. 

To answer these questions, we consider an arbitrary two-terminal mesoscopic 
conductor characterized by the distribution of transmission eigenvalues p{T). The CGF 
for this case can be found with the help of the formula obtained by Levitov et al 
By expanding S{\) [equation ()A.1|) in the Appendix] in Taylor series up to O(A^), we 
obtain the third cumulant in the form \ 



Ci = eGj de ih - fn) {l - 37i[/l(1 - /«) + - h)] + 272(/l - //?)'} , (5.1) 



where G = Gq Jq p{T)TdT is the Landauer conductance and the other constants are 



I The first and the second cumulant of current fluctuations for arbitrary conductor are given by 
equations (|4.5(l and (|4.6(l . respectively, in which G and F are determined by p{T). 




(5.2) 
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with a = Gq/G. Integrating over energy under the assumption max(fcT, eV^) <^ fik, 
k = L,R, we obtain 

Ci = ekTG [p [cothM + u{l - coth^w)] + 2^3^} (5.3) 

as a function of the bias-to-temperature ratio u = eV/{2kT). Asymptotically, 

fe^JF, eV^kT 
^ \ eG[{(3kT + F^eV], eV ^ kT, ^ ' ^ 

where 

F = 1 - 7i = a pT(l - T)dT, 
Jo 

(3 = 6(71 - 72) = 6a /' pT\l - T)dT , (5.5) 

JO 

Fa = 1-371 + 272 = a / pT(l-r)(l-2r)dr. 

Jo 

It is seen that at low biases the third-order cumulant is given by the universal relation 
Gi = e^IF for any distribution of transmission eigenvalues p{T). This means that, at 
low biases, Gi for any conductor is: (i) independent of temperature, (ii) linear in voltage 
with a positive slope, and (iii) proportional to the Fano factor F. At high biases, the 
voltage dependence of C/ is also linear, with the slope determined by the constant F3. 
The slopes in the low-bias and high-bias regimes are different except in the unique case 
of a tunnel junction, for which all the tunneling probabilities are small, T -C 1, and 
F F'i 1, (3 ^ Q. The latter is in agreement with the prediction by Levitov and 
Reznikov [H]. 

The slope of Gj can be negative, that means C/(K) changes sign, if and only if 
F3 < 0. Let us focus our analysis on the sign of the constant F3. We observe that the 
function ip{T) = T(l — T)(l — 2T) obeys ?/'(! — T) = —iIj{T). Hence we can write 

F, = a ["\p{T)-p{l-r)]ij{T)dT, (5.6) 

JO 

where the range of integration is reduced to < T < 1/2. It is clear that the sign 
of this integral depends on the function in the square brackets only, since ip{T) > 0. 
Denoting Ap(T) = p(T) — p(l — T), we observe that Ap{T) is related to the asymmetry 
or skewness of the distribution p(T) around the midpoint T=l/2. If the backscattering 
dominates, we can expect Ap(T) > and F3 > 0. If the forward transmission 
dominates, one can expect Ap(T) < and F3 < 0. This effect is much better pronounced 
in the case of a one-channel conductor jEj, for which p(T) = 6(T — Ti) and the sign of 
F3 depends on the condition whether Ti is greater or smaller than 1/2. 

In principle, one can calculate the coefficient -F3 for any mesoscopic conductor, for 
which the distribution p(T) is known (for different types of the distributions see, e.g., 
the reviews |H]). We just remark that for diffusive metals the distribution function 
is bimodal p(T) ~ l/{T\/\ — T) with clear asymmetry in favor of backscattering. 
One can verify that Ap > and therefore C/ does not change sign. Indeed, direct 
calculation from equations (j5.2|) and (|5.5|) gives, for diffusive metals, 7i=2/3, 72=8/15, 
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F=l/3, /3=4/5, ^3=1/15. This is in agreement with the studies by Nagaev [T^ 
(semiclassical approach) and Gutman and Gefen [21] (nonhnear a-model), who obtained 
the temperature dependent third cumulant for this case. 

For a double-barrier tunnel junction the distribution of eigenvalues is also bimodal: 
p(T) ~ 1/ \jT^i%nax — T) with T^in < T < T^ax [HHj ■ Again, as in the case of diffusive 
metals, there is a skewness of p{T) in favor of backscattering because of two reasons: 
the shape of the distribution and the gap at high transmissions close to T^ax (the gap 
at small transmissions is much smaller [Sn])- The coefficients for this case can be found: 
71 = 1 - F, 72 = (3/2)(l - F)\ p = -3(1 - F)(l - 3F), F3 = 1 - 3F + 3^^, where 
the Fano factor F = {g\ + gji)/{gL + gnY is determined by the conductances gi and 
gji of the barriers. The coefficient F3 in the zero-temperature limit for this case has 
been derived by de Jong jHOl- It is seen that F3 > for any F, including the case of a 
symmetric junction for which the Fano factor takes its maximal value F=l/2. Therefore 
Cj for double-barrier tunnel junctions is always positive. 

For chaotic cavities, the function p(T) is bimodal again [see equation ()A.4|1 in 
the Appendix]; however, in contrast to the previous cases, it is asymmetric in favor 
of the forward scattering due to a gap at low transmissions < T < 7^. Therefore 
for asymmetric cavities Ap < 0. The gap disappears only when the cavity is ideally 
symmetric. In this case p{T) ~ 1/ ^T{1 — T) and Ap=0. Indeed, from equation 
(|4.1ip follows that -F3 = —(1 — 4F) < 0, which is equal to zero only when F=l/4 
for a symmetric cavity. Therefore the negative slope of the volt age- dependent third 
cumulant in asymmetric chaotic cavities is caused by the skewness in the transmissions 
when forward scattering prevails over backward scattering. The voltage independent 
asymptotics occurs only when the function p(T) is symmetric with respect to the 
midpoint T=l/2. 

Summarizing these examples, we conclude that the third-order cumulant of current 
fluctuations is an important quantity characterizing the asymmetry of the distribution 
of transmission eigenvalues p{T) with respect to forward and backward scattering. 
While the average value is described by the conductance G, and the "bimodality" 
(dispersion) is determined by the Fano factor F, the skewness is captured by the factor 
F3, in particular by its sign. Therefore, the measurements of the third-order cumulant 
would give information additional to that obtained from the first- and the second-order 
cumulants. 

6. Summary 

By using the extended Keldysh-Green function technique of the circuit theory, we 
have derived the analytical formula for the cumulant generating function of current 
fluctuations in a chaotic cavity with asymmetric leads at an arbitrary bias-to- 
temperature ratio. The CGF contains a complete information about all-order cumulants, 
and it is expressed in terms of the anticommutator of the matrix Green functions at the 
leads, and therefore can be applied to different types of the contacts (see footnotef). 
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We emphasize the importance of knowledge of the CGF for any stochastic problem: 
it gives not only the cumulant of any order, but also the probability distribution of 
a certain event including "rare" events (extreme- value statistics), which cannot be 
obtained from the knowledge of the first few cumulants. 

For metal contacts with Fermi-Dirac filling factors and fully transparent leads, going 
beyond the previous studies at zero temperature, we have found that the third cumulant 
of current noise Cj changes sign at high biases for any asymmetric cavity except in the 
case when the number of quantum channels in the two leads are equal (symmetric 
cavity). This effect is attributed to the skewness of the distribution of transmission 
eigenvalues with respect to forward/backward scattering. For the asymmetric case, 
Cj is linear in bias at both low and high voltages with positive and negative slopes, 
respectively. For the symmetric cavity, Ci is linear at low voltages and it is constant at 
high voltages with the value proportional to the temperature and the number of quantum 
channels. This value was found to be related to the cascade correlations between the 
cumulants. At low biases (or, equivalently, at high temperatures), the third cumulant 
is proportional to the Fano factor F and independent of temperature. This fact can be 
used in the experiment to measure the Fano factor F from the ratio Ci / e^I at arbitrarily 
small biases. In contrast to the case of tunnel junctions, for which Cj may change sign 
due to effect of the environment, in asymmetric cavities Cj changes sign due to internal 
dynamics. The correspondence with other theories: the RMT theory, the cascaded 
Langevin approach, and the stochastic path integral method has been discussed. 

The first successful measurements of the third cumulant in tunnel junctions jU E] 
leads us to believe that similar measurements can be performed for cavities and other 
mesoscopic systems. 

Appendix A. Correspondence with RMT 

Levitov's formula for the CGF provides a connection between the FCS and the scattering 
properties of a two-terminal mesoscopic conductor 0: 



where p{T) is the distribution of transmission eigenvalues. In the zero-temperature limit 
this formula can be inverted to find the function p(T) from the CGF S{X) by analytical 
continuation in the complex plane [28] 



where A = e^'^ — 1. One can find p{T) for an asymmteric cavity by substituting Eq. 
()2.18|) into ()A.2|) . The shortest way, however, is to use directly the derivative from 
()2.14|1 . one obtains 




(A.l) 




(A.2) 




(A.3) 
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where % = 1 - AaiaR = {Nl - NrY/{Nl + NrY. Now applying equation (TOl) to 
()A.3|) . we arrive at: 

P(T) = ^^^, r.<T<l, (A.4) 

and p{T) = at < T < 7^. This coincides with the distribution obtained by Beenakker 

[H] from the random matrix theory based upon the circular unitary ensemble of chaotic 
cavities. 
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Figure 1. Chaotic cavity with two asymmetric leads supporting Nj^ and quantum 
channels (left) and its circuit theory representation (right) with matrix Green functions 
Gl, Gc, and Gr assigned to each terminal and to a central node. The matrix currents 
II and Ir obey a Kirchoff law in a central node. 




Figure 2. Counting statistics of charge Q transmitted trough an asymmetric cavity 
for several parameters of the asymmetry rj = Nr/N^. The narrowest contact R is kept 
fixed, while the contact L is widened. The charge is normalized to its maximum value 

Qmax = {s'^yt/'KK)Nii determined by the conductance of the narrowest opening R, 
and the probability is normalized according to P{Q) — 1 at the point of the average 
transmitted charge Q — QlQr/{Ql + Qr)- 
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Figure 3. Third cumulant of current fluctuations Cj as a function of bias-to- 
temperature ratio for different degrees of the cavity asymmetry i] = Nfi/N^. With 
increasing bias, C/ changes sign for any rj except 77=1 for the ideally symmetric case, 
for which it saturates. For negative biases Ci{—V) — —Ci{V). 
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Figure 4. Voltage dependence of the third cumulant Ci (blue) for the asymmetric 
cavity with 77 = 1/3. At low and high biases Cj is approximately linear (red) with 
positive and negative slopes, respectively. 



